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An analytic method for deriving the free energy of a three-dimensional Ising- 

like system near the critical point in a homogeneous external field is developed in 

the p^ model approximation. The mathematical description proposed for temper- 

2 \ atures T > Tc {Tc is the phase-transition temperature in the absence of an external 

field) is valid for fields near he, where the scaling variable is of the order of unity 

'^ ■ and power series in this variable are not effective. At the limiting field he, the 

^ , temperature and field effects on the system in the vicinity of the critical point are 

Q ! equivalent. The total free energy is obtained as a function of temperature, field 

and microscopic parameters of the system without using series expansions in the 

scaling variable. In addition to leading terms, the expression for the free energy 

vQ ■ includes the terms determining the temperature and field confluent corrections. 
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i> ■ 1 Introduction 

o 

In spite of great successes in the investigation of three-dimensional {SD) 
k> ! Ising-like systems made by means of various methods (see, for example, [1]), 

H i the study of the effect of an external magnetic field on the critical behaviour 

of the mentioned systems and the calculation of scaling functions are still of 
interest [2] . The Ising model is widely used in the theory of phase transitions 
for the study of the properties of various magnetic and non-magnetic systems 
(ferroelectrics, ferromagnets, binary mixtures, etc.). 

In this article, the behaviour of a 3D Ising-like system near the critical 
point in a homogeneous external field is studied using the collective variables 



(CV) method [3-5]. The main pecuharity of this method is the integration 
of short-wave spin-density oscillation modes, which is generally done with- 
out using perturbation theory. The CV method is similar to the Wilson 
non-perturbative renormalization-group (RG) approach (integration on fast 
modes and construction of an effective theory for slow modes) [6-8]. The 
term collective variables is a common name for a special class of variables 
that are specific for each individual physical system [3,4]. The CV set con- 
tains variables associated with order parameters. Because of this, the phase 
space of CV is most natural for describing a phase transition. For magnetic 
systems, the CV pk are the variables associated with modes of spin-moment 
density oscillations, while the order parameter is related to the variable po, 
in which the subscript "0" corresponds to the peak of the Fourier transform 
of the interaction potential. 

The free energy of a 3D Ising-like system in an external field at tem- 
peratures above Tc is calculated using the non-Gaussian spin-density fluctu- 
ations, namely the sextic measure density. The latter is represented as an 
exponential function of the CV whose argument includes the powers with the 
corresponding coupling constants up to the sixth power of the variable (the 
p^ model). The present publication supplements the earlier works [9-12], in 
which the p^ model was used for calculating the free energy and other ther- 
modynamic functions of the system in the absence of an external field. The 
p^ model provides a better quantitative description of the critical behaviour 
of a 3D Ising-like magnet than the p^ model [11]. 

2 Basic relations 

We consider a 3D Ising-like system on a simple cubic lattice with A^ sites 
and period c in a homogeneous external field h. The Hamiltonian of such a 
system has the form 

H = -l E-^lniVj^i-^E^J- (1) 

Here rji is the distance between particles at sites j and 1, and ctj is the 
operator of the z component of spin at the jth site, having two eigenvalues 
+ 1 and -1. The interaction potential is an exponentially decreasing function 
<l>(rji) = Aexp (— Tji/6), where A is a constant and b is the radius of effective 



interaction. For the Fourier transform of the interaction potential, we use 
the following approximation [3,11,12]: 

' $(0)(1-262A;2), k<B', 
0, B' <k< B. 
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Here B is the boundary of the Brillouin half- zone (B = vr/c), B' = {b^/2) ^, 
$(0) = 8TTA{b/cf. 

The integration over the zeroth, first, second, . . ., nth layers of the CV 
phase space [3-5, 11] leads to the representation of the partition function of 
the system in the form of a product of the partial partition functions Qn of 
individual layers and the integral of the "smoothed" effective measure density 

Z = 2^2(^"+-^)/2Qogi ■ ■■Qn[Q{Pn)f"^'Jwt^'\p) (dp)^"+^ (3) 

The expressions for Qn, Q{Pn) are presented in [9-12], and A^„+i = A^'s^'^("'+^) 
(s is the RG parameter), N' = Nsq'^ {d = 3 is the space dimension), sq = 
B/B' = ny/2b/c. The sextic measure density of the {n + l)th block structure 
Wq (p) has the form 
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where 5„+i = 5's-("+i), rf„+i(A;) 
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thermodynamic temperature, a^ and a^i are the renormalized values 
of the initial coefficients after integration over n + l layers of the phase 



space of CV. The coefficients Oi = s ^tn, dn{0) 
dn{k) = dn{0) +2f3^0)b^k% ai"^ 
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s '^"'Tn [appearing in 
s~^^Wn are connected 

with the coefficients of the (n + l)th layer through the recurrence relations 

(RR) 

Ad+2)/2. 
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whose solutions 
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«;. = «;(°) + ci4%^°^i^2" + C2«;S?(n(°))i/2E3" + c,E2 (6) 

in the region of the critical regime are used for calculating the free energy of 
the system. Here 

The quantity q = g/3$(0) determines the average value of the Fourier trans- 
form of the potential /3$(-B„+i, -B„) = /5$(0) — g/s^" in the nth layer (in this 
article, g = (1 + s~^)/2 corresponds to the arithmetic mean value of fc^ on 
the interval (1/s, 1]). The basic arguments hn and a„ are determined by the 
coefficients of the sextic measure density of the nth block structure. The 
intermediate variables rjn and ^„ are functions of hn and On- The expressions 
for both basic and intermediate arguments as well as the special functions 
appearing in equations ([7]) are the same as in the absence of an external 
field (see [9-12]). The quantities Ei in equations (E]) are the eigenvalues of 
the matrix of the RG linear transformation. We have Ei = s^'^^'^'''^ . Other 
the eigenvalues £^2, E^ and £"4 coincide, respectively, with the eigenvalues 
El, E2 and £3 obtained in the case oi h = 0. The quantities /o, v^o and 
ipo characterizing the fixed-point coordinates (t'-"-* = 0, r^°) = — /o/3$(0), 
ti(o) = y9Q(/5$(0))^, w^^^ = fpoiP^iO))^) as well as the remaining coefficients 
in equations (jHl) are also defined on the basis of expressions corresponding to 
a zero external field. 

3 Free energy of the system at T > Tc as func- 
tion of temperature, field and microscopic 
parameters 

The basic idea of the free-energy calculation on the microscopic level consists 
in the separate inclusion of the contributions from short-wave (the region 



of the critical regime) and long-wave (the region of the limiting Gaussian 
regime) modes of spin-moment density oscillations [3-5]. The contributions 
from short- and long- wave modes to the free energy F = — kT In Z in the 
presence of an external field are calculated in the p^ model approximation 
according to the scheme proposed in [9-12]. Short-wave modes are charac- 
terized by a RG symmetry and are described by the non-Gaussian measure 
density. The calculation of the contribution from long-wave modes is based 
on using the Gaussian measure density as the basis one. Here, we have de- 
veloped a direct method of calculations with the results obtained by taking 
into account the short-wave modes as initial parameters. 

The main peculiar feature of the present calculations lies in using the 
generalized point of exit of the system from the critical regime of order- 
parameter fluctuations. The inclusion of the more complicated expression 
for the exit point (as a function of both the temperature and field variables) 

- - ^ "^ -1 (8) 
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leads to the distinction between formula for the free energy of the system 
and the analogous relation at /i = [10, 11]. The quantity h = h'/ fo is 
determined by the dimensionless field h' = jSh, while the quantity he = fP° 
is a function of the reduced temperature r = {T — Tc)/Tc. Here f = c^ t/ fo, 
Pq = In ii^i / In ii^2 = {d+2)u/2, c^ characterizes the coefficient Ci in solutions 
([6]) of RR, 1/ = In s/lnE2 is the critical exponent of the correlation length. At 
h = 0, Up becomes rrij. = — In f/ In ii^2 ~ 1 (see [5,10,11]). At T = T^ (r = 0), 
the quantity Up coincides with the exit point Uh = — In/i/lnii^i — 1 [13]. The 
limiting value of the field he is obtained by the equality of the exit points 
defined by the temperature and by the field (m^ = n/j). 
Having expression ([8]) for Up, we arrive at the relations 

Ei"-+i = (h' + hi)-'/', tE^^-"^ = He, He = h'J^'Ch' + hlr'^^^^'K 
E^^^' = H3, Hs = {h' + /i2)^i/(2po)^ 
E^^+i = if^, H^ = [h' + /i2)^2/(2po)^ 

S-("^+l) = (/i2 + /,2)l/(rf+2)^ (9) 

where Ai = — lnii^3/lnii^2 and A2 = — lnii^4/lnii^2- In contrast to He, 
the quantity H3 takes on small values with the variation of the field h (see 
figure [I]) . The quantity He at h ^ and near he is close to unity and series 




Figure 1: Behaviour of quantities H^ and H^ with increasing ratio h/hc for 
the RG parameter s = s* = 2.7349 and the reduced temperature r = 10~^. 



expansions in He are not effective here. The value oi s = s* = 2.7349 (see 
figure caption) corresponds to nulhfying the average value of the coefficient 
in the term with the second power in the effective density of measure at the 
fixed point [5, 10, 11]. 

We perform the calculations with the help of equations ([9]), which are 
valid in the general case for the regions of small, intermediate and large field 
values. The inclusion of E^^ (or i^s) leads to the formation of the first 
confluent corrections in the expressions for thermodynamic characteristics of 
the system. The quantity E^'^ (or H^ is responsible for the emergence of 
the second confluent corrections. We disregard the second confluent correc- 
tion in our calculations. This is due to the fact that the contribution from 
the first confluent correction to thermodynamic functions near the critical 
point (r = 0, /i = 0) is more significant than the small contribution from the 
second correction {li^ + /i^ ^ 1, Ai is of the order of 0.5 and A2 > 2, see [9]). 

Collecting the contributions from all regimes of fluctuations at T > Tc 
in the presence of an external field and using the relation for s~("-p+i) from 
equations (jH]), we obtain the following expression for the total free energy of 
the system: 
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The advantage of the method presented in this article is the possibihty of de- 
riving analytic expressions for the free-energy coefficients as functions of the 
microscopic parameters of the system (the lattice constant c and parameters 
of the interaction potential, i.e. the effective radius b of the potential, the 
Fourier transform $(0) of the potential for k = 0). The coefficients of the 
non-analytic component of the free energy depend on H^. The terms pro- 
portional to ifs in equation flTUl) determine the confluent corrections by the 
temperature and fleld. As is seen from the expression for F, the free energy 
of the system aX h = and f = 0, in addition to terms proportional to f^" 
(or h^^^) and h^^^, contains the terms proportional to f^'^+^^ and h^/^+^^/po^ 
respectively. At h ^ and f 7^ 0, the terms of both types are present. 

4 Conclusions 

An analytic method for calculating the total free energy of a 3D Ising-like sys- 
tem near the critical point is developed on the microscopic level in the higher 
non-Gaussian approximation based on the sextic distribution for modes of 
spin-moment density oscillations (the p^ model). The simultaneous effect of 
the temperature and fleld on the behaviour of the system is taken into ac- 
count. An external fleld is introduced in the Hamiltonian of the system from 
the outset. In contrast to previous studies on the basis of the asymmetric 
p"^ model [13-16], the fleld in the initial process of calculating the partition 
function of the system is not included in the Jacobian of transition from the 
set of spin variables to the set of CV. Such an approach leads to the appear- 
ance of the flrst, second, fourth and sixth powers of CV in the expression for 
the partition function and allows us to simplify the mathematical description 
because the odd part is represented only by the linear term. 

The main distinctive feature of the proposed method is the separate in- 
clusion of the contributions to the free energy from the short- and long- wave 
spin-density oscillation modes. The generalized point of exit of the system 
from the critical regime contains both the temperature and fleld variables. 
The form of the temperature and fleld dependences for the free energy of 
the system is determined by solutions of RR near the flxed point. The ex- 
pression for the free energy obtained at temperatures T > Tc without using 
power series in the scaling variable and without any adjustable parameters 
can be employed in the fleld region near he- The limiting fleld he satisfles the 
condition of the equality of sizes of the critical-regime region by the tempera- 



ture and field (the temperature and field effects on the system in the vicinity 
of the critical point are equivalent) [13-16]. Proceeding from the expression 
for the free energy, which involves the leading terms and terms determining 
the temperature and field confluent corrections, we can find other thermo- 
dynamic characteristics (average spin moment, susceptibility, entropy and 
specific heat) by direct differentiation of equation fITU]) with respect to field 
or temperature. 
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